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ABSTRACT 

We present observations from X-ray to mid-infrared wavelengths of the most energetic non-quasar transient ever observed, 
AT20211wx. Our data show a single optical brightening by a factor > 100 to a luminosity of 7 x 10% erg s^! , and a total radiated 
energy of 1.5 x 10°? erg, both greater than any known optical transient. The decline is smooth and exponential and the ultra-violet 
- optical spectral energy distribution resembles a black body with temperature 1.2 x 10^ K. Tentative X-ray detections indicate 
a secondary mode of emission, while a delayed mid-infrared flare points to the presence of dust surrounding the transient. The 
spectra are similar to recently discovered optical flares in known active galactic nuclei but lack some characteristic features. The 
lack of emission for the previous seven years is inconsistent with the short-term, stochastic variability observed in quasars, while 
the extreme luminosity and long timescale of the transient disfavour the disruption of a single solar-mass star. The luminosity 
could be generated by the disruption of a much more massive star, but the likelihood of such an event occurring is small. A 
plausible scenario is the accretion of a giant molecular cloud by a dormant black hole of 105 — 10? solar masses. AT20211wx 
thus represents an extreme extension of the known scenarios of black hole accretion. 
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1 INTRODUCTION 


The accretion of matter onto supermassive black holes (SMBHs) is 
the most efficient known process of extracting energy in the Universe, 
with general relativistic effects combined with complex magneto- 
hydrodynamics and highly non-thermal radiative processes causing 
spectacular electromagnetic phenomena. This accretion can happen 
at a widely varying range of rates: steadily-accreting SMBHs, called 
active galactic nuclei (AGN), are the central engines of galaxies, 
which accrete gas for Myr at a time with a rate that typically varies 
on timescales of seconds to years (e.g. McHardy et al. 2006; Scaringi 
et al. 2015). 

As opposed to typical variability of AGN caused by stochastic ac- 
cretion, tidal disruption events (TDEs) represent a much shorter, yet 
far more violent accretion episode. These phenomena are caused by 
the destruction of a star in the vicinity of a SMBH (or intermediate- 
mass black hole; Angus et al. 2022) due to tidal forces. TDEs 
usually reveal themselves through a single flare observed in opti- 
cal/ultraviolet (UV)/X-rays with a smooth rise and decaying light 
curve. Not all SMBHs are equally likely to produce a TDE: for 
SMBHs with masses Mpy z 108 Mo, the tidal radius of a solar 
mass star lies within the innermost stable circular orbit (ISCO) of the 
black hole, which suppresses or completely prohibits TDEs (the ex- 
act limit depends on stellar type and SMBH spin). Indeed, observed 
or model-inferred masses of TDEs suggest Mgy ~ 5x 10? 2 107 Mo 
and M, ~ 0.6—13M o (Mockler et al. 2019; Ryu et al. 2020), though 
the high end of this stellar mass range is highly uncertain. TDEs dis- 
play broad emission lines, often of hydrogen and/or helium, and 
many show ionised iron and Bowen flourescence lines of nitrogen 
and oxygen (e.g. Arcavi et al. 2014; Leloudas et al. 2019; van Velzen 
et al. 2021), indicative of the presence of a UV-bright accretion disk. 

The clear differences between regular AGN variability and TDEs 
are muddied by the addition of recently-discovered changing look 
AGN (CLAGNSs; LaMassa et al. 2015) and SMBH ‘re-ignitions’ 
(Trakhtenbrot et al. 2019). CLAGNs tend to exhibit changes in their 
spectral properties such as line widths and continuum slopes (Ko- 
mossa et al. 2022), sometimes accompanied by X-ray outbursts and 
a change in optical brightness but critically while maintaining their 
usual variable optical light curve. On the other hand, re-ignitions and 
flares appear as a single highly luminous optical transient (Kankare 
et al. 2017; Gromadzki et al. 2019; Trakhtenbrot et al. 2019; Fred- 
erick et al. 2021). The spectra of these flares, which are often found 
in narrow-line Seyfert 1 (NLSy1) host galaxies, show ubiquitous hy- 
drogen Balmer emission lines as well as many of the same helium 
and Bowen features of TDEs (Frederick et al. 2021), although the 
line profiles (particularly of hydrogen) tend to be narrower. While 
the high luminosity of the flares necessitates a large amount of mate- 
rial being accreted, the exact nature remains a mystery, with various 
sources labelled as possible abnormal TDEs or sudden changes in 
the accretion flow of already-accreting SMBHs. 

Here we present multiwavelength observations of AT20211wx, 
the most luminous of such flares ever detected. AT20211wx, at z = 
0.9945, displays properties spanning TDEs, SMBH flares, and AGN, 
while unlike all of those classes there is no host galaxy visible to deep 
limits in pre-outburst survey data. Our observations (Section 2) span 
nearly three years and from MIR to X-ray wavelengths. We present 
details of our photometric (Section 3) and spectroscopic (Section 
4) modelling, motivating our discussion of the potential scenarios 
of the event (Section 5), after which we conclude by interpreting 
AT20211wx as an accretion of a non-stellar, gaseous object onto a 
large SMBH. 


Throughout the paper, we use 'days' to refer to the observer frame 
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Table 1. Summary of the multiwavelength observations of AT20211wx. 


MJD Epoch Instrument Bands 
(rest frame d) 
Photometry 
58750 to 59937 —271 to +323 ZTF gor 
58113 to 59908 —590 to +309 ATLAS c, 0 
56803 to 59966, | —1247to +114, WISE WLW2 
59923, 59966 +317, +338 Swift UVOT uvw2,uvm2, 
uvwl,u,b,v 
X-ray 
59923, 59966 +317, +338 Swift XRT 0.3 — 10 keV 
Spectroscopy 
59345, 59844 +27, +277 NTT EFOSC2  3685-9315À 
59931 +320 GTC EMIR 0.85 — 1.35 um 
59949 +329 GTC EMIR 1.45 — 2.42 um 


and ‘d’ for the rest frame at z = 0.9945. All photometry is cor- 
rected for Galactic foreground extinction according to Schlafly & 
Finkbeiner (2011) assuming a Fitzpatrick (1999) extinction curve. 
Magnitudes are presented in the AB system (Oke & Gunn 1983). 
Upper limits and uncertainties are at the 1o level. Where necessary 
we assume a spatially flat lambda cold dark matter cosmological 
model with Ho = 70 km s^! Mpc^! and Qu = 0.3. 


2 OBSERVATIONS 


AT2021Iwx was discovered as a transient in Zwicky Transient Facil- 
ity (ZTF; Bellm et al. 2019; Graham et al. 2019) imaging on the 13th 
of April 2021. Forced photometry revealed it to have been brightening 
since at least 16th June 2020 as ZTF20abrbeie. It was independently 
detected by the Asteroid Terrestrial-impact Last Alert System (AT- 
LAS; Tonry et al. 2018; Smith et al. 2020) as ATLAS20bkdj and 
the Panoramic Survey Telescope and Rapid Response System (Pan- 
STARRS; Wainscoat et al. 2016; Chambers et al. 2016) as PS22iin. 
AT20211wx was twice observed by the extended Public ESO Spec- 
troscopic Survey for Transient Objects ‘plus’ (ePESSTO+; Smartt 
et al. 2015) and was classified as an AGN at z = 0.995 (Grayling 
et al. 2022). In this section we describe photometric observations 
of AT20211wx at X-ray, ultraviolet, optical, near- and mid-infrared 
wavelengths, which are summarised in Figure 1. Times f are given in 
the rest frame at z = 0.9945 with respect to maximum light, which 
we take as the date of the brightest r-band observation: MJD 59291. 


2.1 Optical photometry 


The optical light curve is compiled using photometry from publicly 
available ZTF g- and r-band photometry!, which we obtained via 
the ZTF forced photometry service (Masci et al. 2019). In addition, 
we obtained ATLAS c (“cyan", 4200 — 6500 À) and o (“orange", 
5600 — 8200 À) imaging over the course of the light curve. The 
region was also covered infrequently as part of the ongoing Pan- 
STARRS Near Earth Object search which provides grizy coverage. 


1 https://lasair-ztf.Isst.ac.uk/object/ZTF20abrbeie/ 
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Figure 1. a) Light curve of AT20211wx. Epochs of our multiwavelength follow-up observations are marked with dashed lines; b) Pan-STARRS upper limits 
up to 750d (rest frame) before the first detection of AT20211wx; c) Comparison to similar transients: NLSyl accretion events AT2019fdr and AT20219brs 
(Frederick et al. 2021) and the prototype nuclear flare PS1-10adi (Kankare et al. 2017); the most luminous known likely TDE ASASSN-151h (Leloudas et al. 
2016); the most luminous known supernova SN2016aps (Nicholl et al. 2020); a likely jetted TDE AT2019dsg (van Velzen et al. 2021), a typical TDE AT2018dyb 
(Leloudas et al. 2019), and the MIR-brightening TDE PS16dtm (Petrushevska et al. 2023) 


Data from both ATLAS and Pan-STARRS were obtained as forced 
photometry (Smith et al. 2020). 

To search for a host galaxy, we use the stacked grizy images 
from the Pan-STARRS 37 survey (Waters et al. 2020). We use the 
photutils package and perform photometry in a circular aperture 
with a 1” radius centered on the location of AT2021]wx, and use 
the PS1 weight image to estimate the photometric uncertainties. 


22 Mid-infrared photometry 


The location of AT20211wx was observed by the Wide-field Infrared 
Survey Explorer (WISE) spacecraft as part of the NEOWISE re- 
activation mission (Mainzer et al. 2011). We obtained NEOWISE 
photometry in the W1 and W2 bands from the NASA/IPAC infrared 


? https://fallingstar-data.com/forcedphot/ 


science archive?. Observations of the location exist with a six month 
cadence from several years before the flare until November 2021, dur- 
ing its decline. In particular we highlight MJDs 58993 (pre-flare), 
59156 (on the rise), 59357 (around peak), and 59520 (during the 
decline). For each epoch, we use the IRSA WISE/NEOWISE coad- 
der (Masci & Fowler 2009) to combine individual frames. We use a 
sigma-clipped median to estimate the background, which we subtract 
and perform aperture photometry in a 4" radius around the tran- 
sient location, corresponding roughly to the point spread function of 
WISE. Before the optical transient, we find no significant detection; 
after the optical flare begins there is a clear point source in both W1 
and W2. 


? https://irsa.ipac.caltech.edu/Missions/wise.html 
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Figure 2. The rest frame UV to MIR SED of AT20211wx at +338 d. The MIR 
data come from +114 d and are not scaled as we do not predict their evolution. 
Data are corrected for MW reddening, but not for any host extinction. The 
uvw2 point is not included in the black body fit due to likely Lyman-a@ 
absorption. 


2.3 X-ray and ultraviolet observations 


We obtained two epochs of observations with the Neil Gehrels Swift 
Observatory (Swift; Gehrels et al. 2004). The first observation took 
place at MJD 59923.554 (tmax+317 d), and lasted 2670 s. The second 
took place at MJD 59966.483 (tmax + 338 d) and lasted 2605 s. Data 
were observed in photon counting mode with the X-ray Telescope 
(XRT; Burrows et al. 2005), and in the uvw2 (central wavelength, 
1928 À), uvm2 (2246 À), uvw1 (2600 À), u (3465 À), b (4392 À) 
and v (5468 À) filters with the Ultraviolet-Optical Telescope (UVOT; 
Roming et al. 2005), respectively. 

The XRT data were reduced with the tasks xRTPIPELINE and xs- 
ELECT. The source and background events were extracted using a 
circular region of 40" and an annular ring with inner and outer 
radii of 60" and 110”, respectively, both centered at the position 
of the source. The Ancillary Response Files (ARFs) were created 
with the task xRTMKARF and the Response Matrix File (RMF), 
swxpcOto12s6 20130101v014.rmf, was taken from the Calibration 
Data base^. Due to the low count rates, the XRT spectra were grouped 
to have a minimum of 3 counts per bin using the FTOOL GRPPHA. 

In terms of the UVOT data, we used the task uvoriMsUM to sum 
all the exposures when more than one snapshot was included in each 
individual filter data and the task uvorsource to extract magnitudes 
from aperture photometry. A circular region of 5" centered at 
the target position was chosen for the source event and another 
region of 40" located at a nearby position was used to estimate the 
background emission. 


2.4 Optical spectroscopy 


Optical spectra were obtained on MJD 59345 and MJD 59844 
(tmax + 27 d and tmax + 277 d respectively) via the European South- 


4 https://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/swift 
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ern Observatory (ESO) as part of ePESSTO+ using the ESO Faint 
Object Spectrograph and Camera (EFOSC2; Buzzoni et al. 1984) on 
the New Technology Telescope (NTT) at ESO La Silla observatory, 
Chile. The first epoch consisted of a single 600s exposure, while 
the second epoch consisted of two 2700 s exposures. All three spec- 
tra made use of grism13 (3685 — 9315 A, resolution ~ 17 A). These 
spectra were reduced using the PESSTO pipeline (Smartt et al. 2015) 
v3.0.1. 


2.5 Near infra-red spectroscopy 


We observed AT2021Iwx on MJD 59931 and MJD 59949 (tmax + 
320d and tmax + 329d respectively) with the Espectrógrafo Mul- 
tiobjeto Infra-Rojo (EMIR, Balcells et al. 2000) on the 10.8m 
Gran Telescopio Canarias (GTC) at Observatorio del Roque de 
los Muchachos, La Palma, Spain. Spectra were obtained under Di- 
rector’s Discretionary Time proposal GTC05-22BDDT (PI: Müller 
Bravo) in long-slit mode. The first epoch made use of the YJ grism 
(0.85 — 1.35 um, resolution ^ 10 A); the second epoch used the HK 
grism (1.45 — 2.42 um, resolution ~ 20À). 

Both epochs used a 1.0" slit. Image rectification, dark and flat 
fielding, sky subtraction, dither combination, and wavelength cali- 
bration were performed using PyEMIR (Pascual et al. 2010; Cardiel 
et al. 2019) v0.17. Cosmic rays were detected and removed using 
the lacosmic algorithm (van Dokkum 2001). Trace extraction was 
performed manually using python. Telluric features were removed 
by observing and subtracting the spectrum of the telluric standard 
star HIP104599. 


3 PHOTOMETRIC PROPERTIES 


The light curve of AT20211Iwx (Fig. 1a) peaked in the observer-frame 
r-band on MJD 59291, 122 d after and 3 mag brighter than the first 
ZTF detection. The decay tends to a smooth decline of ~ 0.004 mag 
d^, consistent with an exponential (e 9.9057) decay, while the late 
time decline is also consistent with the power-law 13/5. To measure 
the (pseudo-) bolometric light curve we K-correct the observed mag- 
nitudes according to Hogg et al. (2002) by assuming a black body as 
the underlying spectral energy distribution (SED). By correcting the 
observed r-band to the rest-frame u we are able to compare to similar 
events at lower redshift. For those objects for which rest-frame light 
curves are not public, we estimate a K—correction to the u-band by 
choosing the observer frame band that is closest to u when shifted 
to the rest frame. We then add the correction of —log(1 + z) to ac- 
count for the shortening of the wavelength range of the emission 
in the rest frame. The rest-frame UV peak absolute magnitude is 
M ~ -25.7 mag. 


3.1 Comparison to similar events 


We build a comparison sample by selecting the most luminous ex- 
amples of transients from a variety of observational classes, the rest 
frame u-band absolute magnitude light curves of which we show in 
Fig. 1b. Our choice of objects serves two purposes: first, we search 
the literature for the brightest known objects of various classes; sec- 
ond, we show typical examples of TDEs, the most likely candidate 
for the type of event that could produce such extreme luminosity. 
The light curve and luminosity of AT20211wx most closely resem- 
ble a handful of transients detected in known narrow-line Seyfert 1 
(NLSy 1) galaxies presented by Frederick et al. (2021), from which 
we select the two that display the brightest luminosity and longest 
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Figure 3. Rest frame UV and optical spectra of AT20211wx. For clarity, the GTC EMIR +320 d spectrum (dark blue) has been scaled to match the optical NTT 
EFOSC2 spectrum from +275 d. Common emission features are shown in gold while absorption features are marked in purple. Telluric absorption has been 


masked in grey. 


timescales, AT2019brs and AT2019fdr. To these we add the even 
more luminous PS1-10adi (Kankare et al. 2017). NLSy1s are highly 
accreting AGN, and the flares have been attributed to sudden en- 
hancements of accretion rate (either due to a TDE, or as increased gas 
flow) onto already-accreting SMBHs. At peak, AT20211wx is nearly 
2.5 times more luminous than any of those events. AT20211wx is 
several orders of magnitude brighter than any known supernova (SN; 
De Cia et al. 2018; Angus et al. 2019; Nicholl et al. 2020). We show 
the most luminous, SN2016aps (Nicholl et al. 2020), which has been 
interpreted as a possible pulsational pair instability SN (PPISN). 
AT2021lwx is ten times more luminous than the extreme TDE 
ASASSNIS-Ih (Dong et al. 2016; Leloudas et al. 2016) and three 
times brighter (in the optical/UV) than the jetted TDE AT2022cmc 
(e.g. Andreoni et al. 2022; Pasham et al. 2023). We also show three 
further, more regular TDEs (although the class is diverse; Leloudas 
et al. 2019; van Velzen et al. 2021): AT2019dsg which likely in- 
cludes a relativistic jet (van Velzen et al. 2021), the well-studied 
AT2018dyb (e.g. Leloudas et al. 2019), and the MIR-brightening 
PS16dtm (e.g. Petrushevska et al. 2023). AT20211wx is 100 times 
brighter and decays slower than these TDEs. 


32 Spectral energy distribution 


Using the (K-corrected) UV photometry from t = 317d we fit 
the spectral energy distribution (SED) with a black body with a 
Levenberg-Marquardt least squares minimisation. We omit the far- 
thest UV band, uvw2, since there is evidence in the spectrum that a 
strong Lyman-qa absorber is present. The fit is shown in Fig. 2. We 
measure a black body colour temperature of Tc = 1.2 + 0.1 x 10 K 
and radius of Rgg = 1.5+0.3x 1016 cm, while the t = 338 d photom- 
etry is best fit by consistent values of (Tc = 1.3 + 0.08 x 10^ K) and 
(Rgp = 1.340.2x 1016 cm). By integrating the black body spectrum, 
normalised to the r-band peak brightness, we estimate the maximum 
pseudo-bolometric luminosity of the transient to be 7 x 10^? erg sl, 
which is similar to the median luminosity of quasars (Rakshit et al. 
2020). We estimate the total energy radiated by integrating the black 


body with the simplistic assumption that the SED did not evolve 
over the duration of the transient. Up until December 2022 (440 rest- 
frame days from the onset of the event), the transient has radiated 
1.5 x 10°3 erg. Such a high energy release in a relatively short time 
is unprecedented for a transient and is usually only associated with 
continuously accreting SMBHs. Given that the transient becomes 
redder with time, this assumption likely leads to an underestimate of 
the peak and total bolometric luminosity. 


The XRT observations resulted in the tentative detection of X- 
rays in the range 0.3—10 keV. Modelling the X-ray spectrum as a 
power-law, we measure an observed 0.3—10 keV flux of 2.94 + 1.59 x 
jg-^ erg s7! cm7?. Correcting for Galactic neutral hydrogen col- 
umn density 10?! cm"? (HI4PI Collaboration, et al. 2016) this flux 
corresponds to an unabsorbed luminosity of 1.52 x 1045 erg s^! . This 
luminosity is far higher than expected by extrapolating the UV-optical 
black body and indicates the presence of a separate emission region. 
The size of the SMBH inferred from the UV-optical properties pre- 
cludes the X-rays originating as thermal disk emission, while this 
X-ray luminosity is of similar order to the UV-optical luminosity at 
the same epoch. An Loptical/ Lx-ray ~ 1 is consistent with late-time 
observations of TDEs (e.g. Gezari et al. 2017; van Velzen et al. 2019; 
Wevers et al. 2019). 


Emission in the MIR is also not consistent with the UV-optical 
black body. It continues rising after the UV-optical starts to decline, 
indicative either of reprocessing by dust or from a relativistic jet. 
Similar effects have been seen in TDEs (e.g. van Velzen et al. 2016; 
Jiang et al. 2021; Onori et al. 2022; Petrushevska et al. 2023) and 
are interpreted as the echoes of the UV emission from a dusty ‘torus’ 
surrounding the SMBH. From a pseudo-bolometric luminosity of 
7x 1046 erg s^! we follow Jiang et al. (2021) to estimate a sublimation 
radius of ~ 400 light days. This distance should be reflected in a lag 
between the UV-optical and MIR emission, but the lack of MIR 
cadence has so far prohibited laying constraints on the time of peak 
of the MIR. Again assuming a black body, we measure a colour 
temperature of Tc ~ 10? K and a radius of Rpg = 10!9 cm. 
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Figure 4. Rest frame UV and optical spectra of AT20211wx as well as the same objects compared in Fig. 1. Spectra have been shifted and scaled in order 
to provide a qualitative comparison. Telluric features are marked with ©. Spectra have been retrieved from the WISeREP spectral database(Yaron & Gal-Yam 
2012). 
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Figure 5. Fits to the Balmer lines H£ (left) and Ha (right) in the +321 d spectrum. We also show the locations of [O m1] and [N ri]. The best fit [O 111] fluxes are 
consistent with being caused by noise, while the [N 11] lines are marginally significant. 


3.3 Host galaxy limits 


From the aperture photometry outlined in Section 2.1 we find no 
significant detections in any band, leading to magnitude limits g > 
23.3, r > 232, i > 23.1, z > 22.3, y > 21.4. The i-band limit 
corresponds roughly to a rest-frame g-band absolute magnitude limit 
of M z -21.5 mag. Assuming a mass-to-light ratio M/L = 2 results 
inan stellar mass upper limit of 7x 1019 M. That is, an average Milky 
Way-like galaxy would be marginally detected in the imaging. Galaxy 
stellar mass (M+) is correlated with the mass of central supermassive 
black holes. A rough Mgy/M,, fraction of 0.025 per cent (Reines & 
Volonteri 2015) implies a BH of mass 1 x 107 Mos Mpu < 3x 
108 Mo fora galaxy with the mass of our estimated upper limit. The 
PSI catalogue (Flewelling et al. 2020) contains an object located ~ 3” 
northwest of the location of AT20211Iwx (PSO 318.4511+27.4315), 
with r — 19.25. If we assume this object is at z — 0.9945, then 
the physical separation is 30.8 kpc, meaning it is unlikely to be 
associated with the transient. Using the PS1 magnitudes and their 
uncertainties we estimate a photometric redshift of the nearby object 
with EAZY (Brammer et al. 2008) in its default configuration. We find 
a best fit photometric redshift Zphot = 2.05 + 0.1, which translates 
to a restframe UV brightness of —23 mag. This object is either a 
bright high redshift galaxy, quasar, or a Milky Way star, and is likely 
unrelated to AT20211wx. 


4 SPECTRAL PROPERTIES 


The rest frame UV-optical spectra of AT20211wx (Fig. 3) are domi- 
nated by Balmer emission lines superimposed on a blue continuum 
which becomes redder with time. Low-ionisation species common to 
active galactic nuclei (AGN) including Mg 11 and He 1 as well as semi- 
forbidden C mı] and C ri] are present, but various forbidden nebular 
lines such as [O 11] and [Om] are absent. The spectra are compared 
to similar objects in Fig. 4 and show similarities and differences to 
TDEs, SNe and AGN-flares. Particularly noticeable is the absence of 
lines indicative of Bowen flourescence, which appear to be ubiqui- 


tous in AGN flares and are also common in TDEs. Ha and Hf are 
well represented by a two-component Gaussian profile (Fig. 5). We 
fit the two lines both simultaneously (although the fit is dominated 
by Ha which lies in a much higher signal-to-noise ratio region) and 
independently. The broad component of Hf is only marginally de- 
tected. To ensure that the broad HG component is taken into account, 
we fix it to the same width as that found from Ha but allow the 
normalisation to vary freely. The fit results show the lines comprise a 
narrow component with full-width at half maximum (FWHM) = 430 
kms! , and a broad component with FWHM 5400 km gi , the broad 
component being blueshifted by ~ 500 — 800 km s7}, The width and 
blueshift are largely dependent on the estimation of the continuum 
which is affected by the line corresponding to the end of the spectral 
coverage of the instrument. If the blueshift is confirmed, it indicates 
outflowing material. 


To estimate a reddening we assume Case B recombination for 
narrow Balmer lines (Osterbrock 1989), with a temperature of T = 
10, 000 K and a density of ne = 100 cm ^? and a Calzetti et al. (2000) 
extinction law (all standard choices for star-forming galaxies; we 
note that if the narrow Balmer lines are from the transient then Case 
B may not be appropriate), we measure a modest dust reddening 
of E(B — V) = 0.23. The lack of forbidden oxygen lines in the 
spectrum allows us to place a limit on the star-formation rate (SFR) 
of the host galaxy. Correcting the spectrum for Milky Way (MW) and 
host reddening, we estimate an upper limit on the [O n] luminosity of 
5.6x 10? erg s^! . Using the L([O n]) - SFR relation from Kennicutt 
(1998), modified according to Kewley et al. (2004), we place an upper 
limit of SFR < 3.7 Mo yr-!. Assuming our upper limit on stellar 
mass of 7 x 10!° Mo, the SFR - stellar mass relation for average 
star-forming galaxies at z = 1 predicts SFR ~ 26 Mo yr! (Zahid 
et al. 2012). We conclude that the host of AT20211Iwx is not a highly 
star-forming galaxy. 


The rest-frame UV spectra of AT20211wx display unresolved ab- 
sorption features from singly-ionised metal species including Fe n, 
Mg u and possibly doubly-ionised Al ur. The uvw2 luminosity is a 
factor of 10 smaller than the rest of the UV-optical. In the rest frame 
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at z = 0.9945, uvw2 lies bluewards of Lyman-a. We infer that there 
is a signifcant Lyman-a absorber located at the systemic redshift, in- 
dicating a large reservoir of warm gas, often observed in quasar (e.g. 
Viegas 1995; Wolfe et al. 2005; Péroux et al. 2006) and gamma-ray 
burst (e.g. Krühler et al. 2013; Wiseman et al. 2017) spectra. 


5 ORIGIN AS AN EXTREME ACCRETION EVENT 


AT2021lwx is a highly energetic event, arguably the most lumi- 
nous optical transient ever observed. Such high luminosity limits the 
available mechanisms with which to explain the event, which we 
summarise in this section. 


5.1 Tidal disruption event 


In this section we discuss the possibility that AT20211wx is a TDE. 
We begin by estimating some parameters of the system under the 
assumption of basic accretion physics and with many simplifications, 
before discussing the results of two different modelling approaches. 

A TDE can only occur when the tidal radius lies outside the event 
horizon of the black hole. That tidal radius depends upon the density 
of the star being disrupted such that for a given stellar mass and radius 
there is a corresponding upper black hole mass limit, the Hills mass 
(Hills 1975), above which a TDE cannot occur. For a non-rotating 
black hole, the Hills mass for a typical 1M main-sequence star is 
~ 8x 107 Mc. Assuming a corresponding radiative efficiency of 0.1 
(e.g. Marconi et al. 2004; Alexander & Hickox 2012; Nicholl et al. 
2022), the peak luminosity of AT20211wx corresponds to a mass 
accretion rate of ~ 1.2 Mo yr`!. This accretion rate lies at half of the 
Eddington limit for a black hole of mass ~ 108 Mo, which is far larger 
than the typical black hole mass inferred for TDEs (Mockler et al. 
2019; Wevers et al. 2019; van Velzen et al. 2021) and marginally 
above the Hills mass for a 1Mo star. Assuming an accretion rate 
typical of quasars (M / Meqa = 0.1) implies Mpy ~ 10? Mo. This is 
beyond the Hills mass for stars up to 200 Mo (unless the black hole 
is rotating with a large spin, in which case the Hills mass increases; 
Leloudas et al. 2016) and lies close to the mean black hole mass 
for quasars at z = 1 which is 1085 Mo (McLure & Dunlop 2004). 
Furthermore, at an accretion rate of 1 — 2 Mo yr-!, over two solar 
masses must already have been accreted, placing that mass as a lower 
limit on any star to have been disrupted. 

The energetics of the event are compatible with a TDE if we expand 
the range of progenitor stars to far higher masses than have been 
observed before. Fitting the light curve using Modular Open Source 
Fitter for Transients (MOSFIT; Guillochon et al. 2018) using the TDE 
model of Mockler et al. (2019) we estimate Mg = 8.3 x 108 Mo 
and M, = 14.8 Mo, consistent with accretion close to the Eddington 
limit. It has previously been assumed that TDEs of such massive stars 
are extremely rare, not least because the lifetimes of such massive 
stars are of the order 15 Myr (Meynet & Maeder 2002). A star of 
this mass must therefore be born within the very central region of the 
galaxy in order to pass within the tidal radius of the SMBH within 
its lifetime, which we find unlikely. Similarly, the existence of a 
~ 15 Mo star requires strong star formation in the host galaxy which 
is generally inconsistent with our upper limits. Nevertheless, the MW 
Galactic centre hosts a population of young (x 10—100 Myr) massive 
(= 10 Mo) stars, despite there still not being a clear understanding of 
their origin (e.g. Paumard et al. 2006; Lu et al. 2008). Recent work 
has also suggested that there are more tidal disruptions of higher- 
mass stars than previously expected — possible explanations include 
higher SFR in the centres of these galaxies compared to elsewhere 
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in the galaxies or top-heavy initial mass functions in galactic nuclei 
producing more higher-mass stars in these regions (Mockler et al. 
2022). Furthermore, the apparent overabundance of carbon compared 
to oxygen is consistent with some observations of evolved asymptotic 
giant branch stars (e.g. De Beck & Olofsson 2020). Even so, the lower 
limits from this study were in the 1 — 2 Mo range, which is far below 
the mass inferred for AT2021lwx. Furthermore, we note that there 
is a strong degeneracy in MOSHIT between the radiative efficiency 
and the stellar mass of the accreted star (Mockler & Ramirez-Ruiz 
2021). 


It should also be noted that MOSHIT has not been designed to 
model disruptions of massive stars; it currently only includes zero- 
age main sequence stars, while a massive star may well be evolved into 
a more diffuse, giant stage. Indeed, the disruption of evolved stars, 
or those with diffuse envelopes, is a plausible mechanism for a TDE 
(e.g. MacLeod et al. 2012; Law-Smith et al. 2017), and while uncer- 
tainty in the mass-radius relation from stellar evolution is included 
in MOSFiT’s systematic error measurements, an extreme example 
could push these errors to the maximum. 


Instead of MOSHIT, which relies on simulations of disruption and 
assumes a set of standard relations to describe the transfer of energy 
from the disrupted material, through the accretion disk, into radiation, 
we attempt to model the light curve using specific disk models that 
assume an initial disk condition. Following the work of Mummery 
& Balbus (2020) we are able to explain the peak luminosity and 
late-time UV-optical SED of AT2021Iwx with a 2 Mo star being 
disrupted by a higher mass, maximally spinning SMBH of Mgy > 1x 
10° Mo. These results are reliant on acompact initial disc (around the 
innermost stable circular orbit; ISCO), and a short viscous timescale 
on the order of the orbital timescale — usually disc timescales are 
expected to be several orders of magnitude larger than the orbital 
timescale. We find that including a much larger disc mass (i.e. mass 
of the disrupted star) drives the temperature to be inconsistent with 
observations (although this could be rectified if the observations are 
influenced by significant host galaxy reddening). 


Neither of the above models accounts for the possibility of super- 
Eddington accretion, which is likely to occur in TDEs (Dai et al. 
2018), particularly in the early phases of the event. In the super- 
Eddington regime, the disk is likely to inflate into a geometrically 
thick disk at which point many of the assumptions used in those 
models break down. For example, there will be significant mass loss 
from the disk due to winds, viewing angle effects, and the likely 
launching of relativistic jets. 


The spectra of AT20211wx do not resemble typical TDEs (e.g. 
Leloudas et al. 2019; van Velzen et al. 2021; Charalampopoulos et al. 
2022). Around half of TDEs appear to show He n, Fe m and/or Bowen 
features, and while it is fairly common for TDEs to exhibit only 
Balmer features, such features tend not to include such strong narrow 
components as AT2021]wx (although see PS16dtm; Petrushevska 
et al. 2023 and AT2019dsg; van Velzen et al. 2021). If the transient 
is caused by the disruption of a star, then these spectral features 
must originate from a reprocessing region not present in typical 
TDEs. Instead, the narrow components resemble more closely the 
low-velocity, shocked circumstellar material that gives rise to the 
narrow lines of type IIn SNe. 


Overall, while the light curve, SED and spectra can be explained 
by a TDE, the nature of such a disruption is highly unconstrained 
and significant advances in the modelling of both the disruption of 
massive stars, and the physics of the subsequent accretion, must be 
made before strong conclusions are drawn. 


5.2 Turn on of an AGN via sudden accretion of gas 


An alternative to the disruption of a star is the sudden, and isolated, 
accretion of a large amount of more diffuse material. This is distinct 
from the classical CLAGN (e.g. LaMassa et al. 2015; Frederick 
et al. 2019; Graham et al. 2020) which tend to show a change 
of spectral properties. The turn-on scenario has been proposed 
for the smooth flaring events in known AGN such as AT2017bgt 
(Trakhtenbrot et al. 2019) and AT2019brs (Frederick et al. 2021). 
Unlike those transients, there is no evidence for a pre-existing AGN 
in AT20211wx, although a low-luminosity AGN is not ruled out due 
to the larger cosmological distance. 


A handful of optical nuclear flares that do not adhere to typical 
TDE, CLAGN, or SN characteristics have been accumulated over 
recent years (Kankare et al. 2017; Trakhtenbrot et al. 2019; Frederick 
et al. 2021). Optically, AT20211wx is more luminous than all of these 
events, while the light curve shape is similar to those interpreted 
as unusual accretion events. Spectroscopically, the two-component 
Ha profile is similar to the nuclear transient PS1-10adi, with the 
broad component showing a similar velocity to that transient and 
AT2017bgt, a long-lasting optical flare in a known AGN. The width 
of the broad component in AT20211wx is larger than all other similar 
nuclear flares in Frederick et al. (2021). Unlike all transients in the 
comparison sample, the spectrum of AT20211wx shows no Fer or 
features excited by Bowen fluorescence (e.g. N rr). 

The large widths of the broad lines imply emission from a region 
with a high-velocity dispersion, while the existence of narrow lines 
implies a slow-moving component. Such a scenario resembles a tra- 
ditional AGN according to the unification picture (Antonucci 1993), 
where broad lines emanate from the eponymous broad-line region 
(BLR), comprising clouds of gas that follow virialised orbits around 
a supermassive black hole. If the observed broad Balmer emission 
originates from an illuminated virialised BLR, then the radius of 
this region is rgpg ~ 5 X 1016 — 5 x 10! cm for black holes of 
108 Mo- 10° Mo respectively. This inferred BLR size is similar in 
scale to the inferred black body radius Rgp, although it is likely that 
the assumption of a spherical photosphere does not hold if the UV 
emission originates from an accretion disk. 

The narrow line region (NLR) in AGN comprises gas at much 
larger radii, with light-travel times from black hole to NLR of 
102 — 10? yr, too far to have been ionised by the current flare of 
AT20211wx. If the narrow Balmer lines in the AT20211Iwx spectra 
correspond to a narrow line region in its host galaxy, they must have 
been ionised by a previously UV-bright source (i.e. an AGN's ac- 
cretion disk). However, the lack of narrow nebular lines other than 
from hydrogen and carbon renders the NLR unlikely as the source of 
these features. A second possible explanation for the lack of [O ri], 
[O ur] and narrow Mg ir is the Baldwin effect (Baldwin et al. 1977), 
an observed anti-correlation between AGN luminosity and emission 
line equivalent widths. Indeed, according to the relationship between 
B-band absolute magnitude and [On] equivalent width in Croom 
et al. (2002) we may not expect to see any [O rr] emission at all. [O r1] 
is found to be lacking in around 5 per cent of all AGN (Schawinski 
et al. 2015). 


5.3 Accretion of a gas cloud by a dormant black hole 


The lack of evidence for an AGN coupled with the difficulties of 
explaining the event with a TDE, we explore an alternative origin 
of AT20211wx. It is plausible that a large amount of material in the 
form of a molecular cloud was disrupted (see e.g., Wang et al. 2017) 
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and accreted by the black hole. Assuming a TDE-like scenario, the 
fuelling rate is determined by the fallback rate of material onto the 
black hole after the disruption. To achieve the observed rise time, 
the cloud would need to be compact enough to provide significant 
fuel at early times, since a diffuse cloud would likely result in a 
slower fallback rate. Detailed modelling of this scenario is deferred 
to further work. 


5.4 Lensed superluminous supernova 


The smooth and slow rise, long decline, and colour evolution 
of AT2021lwx resemble superluminous SNe (SLSNe), which are 
thought either to be powered by magnetars (e.g. Woosley 2010; Kasen 
& Bildsten 2010) or interaction of ejecta with dense circumstellar 
material (CSM). The latter scenario is particularly likely in the case 
of SLSNe-IIn (e.g. Smith et al. 2007; Benetti et al. 2014; Nicholl etal. 
2020) which display narrow Balmer features as seen in AT20211wx, 
raising the possibility of a SN origin for the transient. A number 
of observed properties deter us from this conclusion. Primarily, the 
peak luminosity requires an extremely unlikely mass of CSM, of the 
order 250 Mo following Chevalier & Irwin (2011), larger even than 
the zero-age main sequence mass of the largest expected SN progen- 
itors (e.g. Kasen et al. 2011). The required CSM mass is reduced if 
the lightcurve is strongly lensed, and thus magnified, by a foreground 
source. Given the lack of evidence for any galaxy foreground galaxy 
in deep imaging, nor any system at a lower redshift in the spectra, we 
find the lensing explanation unlikely. 


6 CONCLUSIONS 


AT20211wx is an extraordinary event that does not fit into any com- 
mon class of transient. With a total radiated energy > 1053 erg, it is 
one of the most luminous transients ever discovered. By collecting 
and analysing multiwavelength photometry and spectroscopy of the 
transient, we conclude that: 


1) the emission is dominated by a black body with a moderate 
temperate of 12,000 K and a large radius of 1016 cm which cools as 
the transient fades 

ii) there are two components to the material, one with a large ve- 
locity dispersion (^ 5400 km s^!) and another slow-moving compo- 
nent (~ 430 km s7!), the fast-moving component potentially forming 
an outflow 

iii) the luminosity is likely powered by the accretion of a large 
amount of gas onto a supermassive black hole with mass > 108 Mo 

iv) while there is no evidence of AGN activity, there is sign of 
possibly a hot corona or jet as well as the presence of large amounts 
of dust 

v) a tidal disruption of a massive star is unlikely due to the small 
chance of such an event observing 

vi) the spectral and photometric features of the transient suggest 
the sudden accretion of a large amount of gas, potentially a giant 
molecular cloud. 


Further follow-up and modelling of AT20211wx is necessary to reveal 
more about the scenario that caused the flare, and the community is 
strongly encouraged to search for similar events in both the future as 
well as in archival data. 
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